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ABSTRACT 

Stereology allows shifting from the 3D distribution of the volumes of Poisso- 
nian Voronoi Diagrams to their 2D cross-sections. The basic assumption is 
that the 3D statistics of the volumes of the voids in the local Universe has 
a distribution function of the gamma-type. The standard rule of conversion 
from 3D volumes to 2D circles, adopting the standard rules of stereology, pro- 
duces a new probability density function of the radii which contains the Meijer 
G-function. A non-Poissonian distribution of volumes is also considered. The 
distribution of the 3D radii of the Sloan Digital Sky Survey Data Release 7 is 
best fitted by a non-Poissonian distribution in volumes as given by the Kiang 
function with argument of about two. 
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1. INTRODUCTION 

The astronomical analysis of the cellular nature of the large scale structure 
of our universe started with the second CFA2 redshift Survey which produced 
slices showing that the spatial distribution of galaxies is not random but 
is organized in filam e nts w hich represent the 2D projection of 3D bubbles, 
see iGeller fc Huchral (|1989) . The organization of astronomical observations 



contin ued with the 2dF Galaxy Redshift Survey (2dFGRS ). seelColless et al 
( 200lh. and with the Sl oan Digital Sky Survey (SDSS), see I York et all (boooh 



Abazaiian et al (2009). These catalogs of slices allow the determination of 
the size of the voids as approximated by circles of a given radius. A visual 
inspection of these slices allows a rough evaluation of the largest void, which 
turns to be ps 34/h Mpc. A refined statistics requires a dig ital version of th e 
radii as given by the catalog of cosmic voids of SDSS R7, see lPan et all ([201 if ). 

A possible approach to the statistics of the se voids i s given by th e Voronoi 
tessellation, after the two historical papers b vlVoronoil dl907i Il908h . 

Following the nomenclature introduced bv lOkabe et all <| 1992T ) . we call the 
intersection between a plane and the Poissonian Voronoi tessellation (PVT) 
V p (2,3). We briefly re call that the first application of the PVT to astrophysics 
is due to iKiangj (|1966l ) . The applications of Voronoi Diagrams to the galaxies 
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started with llcke fc van de Wevgaert I (|l987l) , where a sequential clustering 
process was adopted in order to insert the in i tial seeds, and continued with 



( 19911); Ivan de Wevgaert 



van de Wevgaert fc Ickd ( 19891): IPierrd (119901): Barrow fe Coled (119 90): Cole; 



Jl99ll): Ivan de Wevgaertl Jl99ll): IZaninettil (199ll): 
Ikeuchi fc Turner Jl99ll ) :ISubba Rao fc Szalavl(ll992l):|van de Wevgaert! (1 1994) : 



Goldwirth etaD f|l995h ; Ivan de Wevgaertl (fcooa . 120031 ): IZaninettil |2006l ). An 
u pdated review of 3D Voronoi Diagrams applied to cosmology can be found 



van de Wevgaertl (|2002Ll2003l ). The 3D PVT can als o be applied to identify 



Marinoni et all (l2002l) : lMelnvk et al.l(|2006l) 
Elviv et al.l |2009l ). 



group s of galaxies in the structure of a super-cl u ster, see Ebcling & Wicdcn mannl 
1 1993); Bcrnardeau & van de Wevgaertl Jl996l): ISchaap fc van de Wevgaertl jgOOOl) : 



van de Wevgaert fc Schaad(|2009l ): 



A different approach to the intersections between bubbles and a plane 
is given by stereology, which is the science of the geometrical relationships 
between structures that exists in three dimensions (3D) and their images, 
which are fundamentally two-dimensional (2D). The absence of a probability 
density function (PDF) for the main parameters of the PVT area in 2D and 
the volume in 3D has not allowed the development of a PDF in radii of the 
V^, (2, 3) problem. Th e publication wi t h a re lative test of a new PDF for the cell 
of PVT as given by iFerenc fc Nedal (|2007l) , allows a simple parametrization 
of the cell. The integral connected with the V p (2,3) problem can now be 
expressed in analytical terms rather than numerical. The previous comments 
can also be rewritten in the form of some key questions. 

• Is it possible to derive the probability density function for the radii of 
2D sections in the Poissonian case? 

• Is it possible to obtain an analytical expression for the survival function, 
see Eq. (|2"5|) . of the radii of 2D sections in the Poissonian case? 

• Is it possible to derive analytical results for the radii of 2D sections in 
the case of non-Poissonian seeds or volumes? 

• Can we apply such obtained analytical results to the catalog of cosmic 
voids as given, for example, by the SDSS R7? 

In this paper we analyze in Section [2] the two main PDFs adopted in 
order to model t he cel ls of PVT which are the old but still wi dely used Kiang 
functi on iKiangl (|1966f ) and the recent Ferenc-Neda function IFerenc fc Neda 
(|2007n . Section[3]reviews the probability of a plane intersecting a given sphere 
, the stereological approach, and then insert in the fundamental integral of the 
stereology the cell's radius of the new PDF. Section [4] contains the observed 
statistics of 1054 cosmic voids , a theoretical comparison with the radii of 
PVT and a comparison of the observed survival function of 2dFGRS with our 
survival function as given by the stereology. An example of non Poissonian 
Voronoi Tessellation (NPVT) statistics at the light of the Kiang function is 
given in Section [5j 
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2. THE DISTRIBUTIONS ADOPTED FOR PVT 

We briefly review the PDFs which regulate the main parameters of PVTs: 
area in 2D, and volume in 3D. 



2.1. The Kiang function 



The gamma variate H(x; c) (|Kiang| (| 1966T ) ) is 



H(x;c) 



r( c ) 



{cx) c 1 exp(— cx), 



(1) 



where < x < oo, c >0, and T is the gamma function. The Kiang PDF has 
a mean of 



and variance 



fi = l, 

c 



(2) 



(3) 



In the case of a ID PVT, c = 2 is an exact analytical result and conversely c 
is suppose d to be 4 or 6 for 2D or 3D PVTs, respectively, the so called Kiang 
conjecture iKiang (|1966l ). 



2.2. Ferenc-Neda function 



A new PDF has been recently introduced. iFerenc fc Nedal (|2007h . in order 
to model the normalized area/volume in 2D/3D PVT 



FN(x; d) =C x x~ exp (-(3d + l)x/2), 
where C is a constant, 

V2V3d+l 



C 



2 2 3 / 2d (3d+ iy 3 / 2d r (3/2 d+ 1/2) ' 



(4) 



(5) 



and d(d = 1, 2, 3) is the dimension of the space under consideration. We will 
call this function the Ferenc-Neda PDF; it has a mean of 



H = 1, 



and variance 



3d+ 1 



(0) 
(7) 



The F erenc-Neda PDF can be obtained from the Kiang function (jKians 
(|1966l )) by the transformation 



3d+ 1 



(8) 
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TABLE 1 

THE VALUES OF x 2 FOR THE CELL NORMALIZED 
AREA-DISTRIBUTION FUNCTION IN 2D AND THE CELL 
NORMALIZED VOLUME-DISTRIBUTION FUNCTION IN 3D ; HERE Tj 
ARE THE THEORETICAL FREQUENCIES AND O/ ARE THE 
SAMPLE FREQUENCIES. HERE WE HAVE 25 087 POISSONIAN 
SEEDS IN 2D, 21 378 POISSONIAN SEEDS IN 3D, AND 40 INTERVALS 

IN THE HISTOGRAM. 



dimension 


PDF 


parameters 


x 2 


2D 


H{x;c)(Eq.^) 


c = 3.55 


83.48 


2D 


f(x;d)(Eq.m 


d = 2 


71.83 


3D 


H(x;c) {Eq.^) 


c= 5.53 


93.86 


3D 


f(x;d) (EqM) 


d = 3 


134.15 



2.3. Numerical results 



In the following, we will model the PVT in which the seeds are computed 
through a random process. The x 2 is computed according to the formula 



N 

E 



Ti 



(9) 



where N is the number of bins, Ti is the theoretical value, and Oi is the exper- 
imental value. A first test of the PDFs presented in the previous section can 
be made by analyzing the Voronoi cell normalized area-distribution in 2D and 
no rmalized volume-distri bution in 3D, see Table[TJ In this comparison f(x; d) 
bv lFerenc fc Nedal ([20071 ). the number of free parameters is zero be c ause d = 2 
or d — 3 fixes the distribution. In the case of H(x; c) by iKianeJ (|l966h . wo 
have one free parameter which is fixed by the sample. 



3. STEREOLOGY 



We first briefly review how a PDF f{x) changes to g(y) when a new variable 
y(x) is introduced. We limit ourselves to the case in which y{x) is a one-to-one 
transformation. The rule for transforming a PDF is 



g{y) 



(10) 



Analytical results have shown that sections through D-dimen sional Voronoi 

Mollerl dl989l 



tessellations are not themselves D-l Voronoi tessellations, see 



1994 


); 


Chiu et al. 


(1996) 
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of a plane intersecting a given sphere is proportional to the sphere's radius, 
R. Cross-sections of radius r may be obtained from any sphere with a radius 
greater than or equal to r. We may now write a general expression for the 
probability of obtaining a cross-section of radius r from the whole distribution 
(which is denoted F(R)): 



which is formula (A7) in iBlower et al.l (|2002h . That is to say, f(r) is the 



probability of finding a bubble of radius R, multiplied by the probability of 
intersecting this bubble, multiplied by the probability of obtaining a slice of 
radius r from this bubble, integrated over the range of R > r. A first example 
is given by the so called monodisperse bubble size distribution (BSD) which 
are bubbles of constant radius R and therefore 

F(R) = ^, (12) 

which is defined in the interval [0, R] and 

= 7W=Wr- (13) 



which is defined in the interval [0, R], see Eq. (A4) in I Blower et al.l (|2002h ■ 
The average value of the radius of the 2D-slices is 

r = l/4i?vr, (14) 

the variance is 

a 1 = 2/3 R 2 - 1/16 i?V, (15) 

and finally, 

Skewness = —1.151, Kurtosis = 0.493. (16) 

3.1. PVT stereology 

In order to find our F(R), we now analyze the distribution in effective 
radius R of the 3D PVT. We assume that the volume of each cell, v, is 

v = ^ttR 3 . (17) 

In the following, we derive the PDF for the radius and related quantities 
relative to the Ferenc-Neda function. The PDF as a function of the radius 
according to the rule of change of variables ([T0|) . is obtained from (|4|) on 
inserting d — 3: 

F(R) = ^^e-f^. (18) 
v ; 243 y ' 

The average radius is 

R = 0.6065, (19) 



G 
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TABLE 2 

THE PARAMETERS OF f(r), EQ. (gg, RELATIVE TO THE PVT CASE. 



Parameter 


value 


Mean 


0.4874 


Variance 


0.02475 


Mode 


0.553 


Skewness 


-.5229 


Kurtosis 


-.1115 



and the variance is 



a 2 (R) = 0.00853. 



(20) 



The introduction of the scale factor, b, with the new variable R — R' /b trans- 
forms Eq. dT5]) into 

20 tv Ft' 3 



F(R',b) = 



400000 7r 5 i?' e 



(21) 



243 6 15 

We now h ave F(R) as given by Eg. (ITO1) and the fundamental integral (TIT]) 
as derived in iFerraro fc Zaninettil ([20111 ) , is 

f(r) = 2/3 K Vsm^rGti rr 2 , 

< r < 1, 



5/6,1/6,1/2 
7/3,2/3, 1/3,0, 



(22) 



where K is a constant, 



K= 1.6485, 



(23) 



Meiier ( 


1936, 


1941); 


Olver et al. 



(j2010h . Details on the real or complex parameters of the Meijer G- function are 
given in the Appendix, [A] Table [5] shows the average value, variance, mode, 
skewness, and kurtosis of the already derived f(r). Asymptotic series are 



f(r) 
when 



2.7855 r 
r ->■ 0, 



and 



f(r) - -0.006 (r - 1 
The distribution function (DF) is 

1^35/6^2/3^,2/ mil 



0.136 (r- 

when r 



(24) 



(25) 



1. 



90 



7r 2 r 6 



1,7/6,1/2,5/6 
8/3,1,2/3,1/3,^,0 



DF[r) = 
1 



7T 

< r < 1. 



(26) 
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TABLE 3 

PARAMETERS OF f(x,b), EQ. (gSJ, RELATIVE TO THE PVT CASE. 



Parameter 


6 = 2.051 


b = 34 


Mean 


1. 


W.57Mpc 


Variance 


0.104 


28.62Mpc 2 


Mode 


1.134 


l8.80Mpc 



The already defined PDF is denned in the interval < r < 1. In order to 
make a comparison with a normalized sample which has a unitarian mean or an 
astronomical sample which has the mean expressed in Mpc, a transformation 
of scale should be introduced. The change of variable is r = x/b and the 
resulting PDF is 



f(x,b) 



2 If VTTi 3/- ,04,1 /100 7T 2 X 6 5/6,1/6,1/2 \ As 



< r < b. 



(27) 



As an example, Table [3] shows the statistical parameters for two different 
values of b. Skewness and kurtosis do not change with a transformation of 
scale. 

We briefly recall that a PDF f(x) is the first derivative of a distribution 
function (DF) F(x) with respect to x. When the DF is unknown but the PDF 
known, we have 

(28) 



F(x) = I f(x)dx. 
lo 



The survival function (SF) S(x) is 

Six) = 1-F(x), 



(29) 



and represents the probability that the variate takes a value greater than x. 
The SF with the scaling parameter b is 



1 - 0.01831 3 5 / 6 10 2 / 3 G^^^ 2 



1,7/6,1/2,5/6 
8/3,1,2/3,1/3,^,0 



SF(x,b) = 
1 



7T 

< r < b. 



(30) 



A first application can be a comparison between the real distribution of 
radii of 1^,(2,3), see Fig. [TJ and the already obtained rescaled PDF f(x,b). 
The fit with the rescaled f(x,b) is shown in Fig. [5] and Table 0] shows the x 2 
of three different fitting functions. 




Fig. 1. PVT diagram V P (2, 3) when 789 2D cells generated by 15 000 3D seeds are 
considered. 



TABLE 4 

THE VALUES OF x 2 FOR THE CELL NORMALIZED 
AREA-DISTRIBUTION OF V p (2,3). THE NUMBER OF 2D CELLS IS 
789, THE 3D SEEDS ARE 15 000 AND THE NUMBER OF BINS IN THE 

HISTOGRAM IS 30. 



PDF parameters x 2 

H(x;c)(Eq. @J)) c = 5.8 250.8 
f(x;d)(Eq. @J>) d= 3.53 250.8 



f{x,b)(Eq. p8}) 6 = 2.0514 127 




Fig. 2. Histogram (step-diagram) of PVT V p (2,3) when 789 2D cells, generated 
by 15 000 3D seeds, are considered. The superposition of the f(x,b), Eq. (f28|) . is 
displayed. 



The PDF Ja of the ar eas of V p (2, 3) can be obtai ned from /(r) by means 
of the transformation, see iFerraro fc Zaninettil (|201ll) , 




that is, 



Ia{A) = /(r) 



/ A (A)=0.549^3^iOG^f^^ 



5/6,1/6,1/2 
7/3,2/3, 1/3,0, 



-2/3 



(31) 



(32) 



The already derived Ja(A) has average value, variance, mode, skewness 
and kurtosis as shown in Table O 

Since, for r close to 0, f(r) ~ r from Eq. (|32|) it follows that /a(0) 7^ 0, in 
particular /a(0) = 0.443 and Fig. [3] shows the graph of /a- 

The previous figure shows that sections through 3-dimensional Voronoi 
tessellations are not themselves 2-dimensional Voronoi tessellations because 
Ja (0) has a finite value rather than as does the 2D area distribution; this fact 
can b e considered a numerical demonstration in agreement with IChiu et al 
(jl996l ). The distribution function Fa is given by 



F A 



0.018 3 5/6 10 2/3 G4;g 



100 A 3 
~9~ ~ 



1,7/6,1/2,5/6 
8/3,1,2/3,1/3,^,0 



(33) 



Consider a three-dimensional Poisson Voronoi diagram and suppose it inter- 
sects a randomly oriented plane 7: the resulting cross sections are polygons. 
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TABLE 5 

PARAMETERS OF f A (A), EQ. {32]), RELATIVE TO THE PVT CASE. 



Parameter 


value 


Mean 


0.824 


Variance 


0.204 


Mode 


0.858 


Skewness 


0.278 


Kurtosis 


-0.337 



i — i i i i i i i — i — i — i — i — i— i— i— i— i— i— i — r 

S N 




J i i i i I i i i i I . i i i I i i i i I i i i i I i i i i L 

0.5 1 1.5 2 3 

A 



Fig. 3. The PDF f A , Eq. (J32J) , as a function of A (full line) and FN(x;d), Eq. (0J, 
when d=2 (dotted line). 
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Fig. 4. Comparison between data (empty circles) and theoretical curve (continuous 
line) of the distribution of areas of the planar cross sections. 

A comparison between Fa and the area of the irregular polygons is shown 
in Fig.|U In this case the number of seeds is 300000 and we processed 100168 
irregular polygons obtained by adding together results of cuts by 41 triples 
of mutually perpendicular planes. The maximum distance between the two 
curves is d max = 0.039. 

As concerns the linear dimension, in our approximation the two-dimensional 
cells were considered circles and thus, for consistency, the radius r of an ir- 
regular polygon was defined as 

A\ 1/2 

that is, r is the radius of a circle with the same area, A, as the polygon. 
The assumption of sphericity can be considered an axiom of the theory here 
presented, but for a more realistic situation the stereological results will be 
far more complex. 

4. STATISTICS OF THE VOIDS 

This section first processes 1024 observed cosmic voids and then derives 
the same results from the stereological point of view. 

4.1. Observed statistics 

The distribution of the effective radius and the radius of the maximal 
enclosed sphere between galaxies of t he Sloan D i gital S ky Survey Data Release 



7 (SDSS DRY) has been reported in iPan et al.l (|201ll) . This catalog contains 



1054 voids: Table [B] shows the basic statistical parameters of the effective 
radius, and Table [3 the radius of the maximal enclosed sphere. 
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TABLE 6 

THE STATISTICAL PARAMETERS OF THE EFFECTIVE RADIUS IN 

SDSS DR7. 



parameter 
elements 
mean 
variance 

standard deviation 

skewness 

kurtosis 

maximum value 
minimum value 



value 
1024 
18.23ft- 1 Mpc 
23.32ft- 2 Mpc 2 
4.82ft- 1 Mpc 
0.51 
0.038 
34.12ft- 1 Mpc 
9.9ft- 1 Mpc 



TABLE 7 

THE STATISTICAL PARAMETERS OF THE RADIUS OF THE 
MAXIMAL ENCLOSED SPHERE IN SDSS DR7. 



parameter value 
elements 1054 

mean 12.95ft" 1 Mpc 

variance 6.99" 2 Mpc 2 

standard deviation 2. 64ft- 1 Mpc 
skewness 1.47 
kurtosis 2.401 

maximum value 25.69ft" 1 Mpc 

minimum value 10ft" 1 Mpc 
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R (Mpc/h) 



Fig. 5. Histogram (step-diagram) of the effective radius in SDSS DR7 with a 
superposition of the PDF in radius of the PVT spheres, F(R, b), as represented by 
Eq. JUJ). The number of bins is 30, and b = 30.05 Mpc. 



4.2. PVT statistics 

Fig.[5]shows a superposition of the effective radius of the voids in the SDSS 
DR7 with a the curve of the theoretical PDF in the radii, F(R,b), as given 
by Eq. ([2~Tj) . Table [8] shows the theoretical statistical parameters. 

Table [9] shows the values of x 2 for the main PDFs here considered. The 
statistics of the voids can also be visual i zed through the SF, see an application 
to the 2dFGRS as given bv lPatiri et al.l (|2006l ): lvon Benda-Beckmann fc Muller 
d2008h . 

The statistics of the voids between galaxies have been also analysed in 
von Benda-Beckmann fc Muller ( 20081 ) with the following self-similar SF in 
the following, Sss, 

1 R \"1 1 R \ p 2 

S ss = e y**) {** x ) , (35) 

where A is the mean separation between galaxies, s\ and S2 are two length fac- 
tors, and pi and pi two powers. A final comp arison between the four samples 
of void size statistics as represented in Fig. 4 of lvon Benda-Beckmann fc Muller 
(|2008r j and our survival function of the radius for V p (2, 3) as given by Eq. ([3"T|) 
is shown in Fig. |6] 

More details as well the PDF of the self-similar distribution can be found 
Zaninettil |201fj|) . 



5. NPVT STATISTICS 



An example of non NPVT is represented by a distribution in volume which 
follows a Kiang function as given by Eq. ([T]). The case of PVT volumes 
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TABLE 8 

THE STATISTICAL PARAMETERS OF THE THEORETICAL RADIUS 
OF THE PVT SPHERES AS REPRESENTED BY EQ. (J2TJ) WHEN 

b=30.05 Mpc. 



parameter value 

mean 18.23ft" 1 Mpc 

variance 7.70ft, -2 Mpc 2 
standard deviation 2.77ft -1 Mpc 

mode 18.49ft -1 Mpc 

skewness 0.0142 

kurtosis —0.0514 



TABLE 9 

VALUES OF x 2 FOR THE EFFECTIVE RADIUS IN SDSS DR7 FOR 
DIFFERENT DISTRIBUTIONS WHEN THE NUMBER OF BINS IS 30. 
IN THIS COMPARISON, THE AVERAGED VALUE OF THE 
ASTRONOMICAL RADII IS ONE. 



PDF 


parameters 


X 2 


H(x;c),(Eq.m 


c= 14.24 


53 


f(x;d),(EqM) 


d = 9.l 


53 


f(x,b),(Eq.m) 


b = 2.051 


182 


F(R,b),(Eq.m) 


b = 1.648 


407 


F K (R,b,c),(Eq.§^ 


6 = 31.33 c= 1.76 


66.121 
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R (h Mpc) 



10 12 



Fig. 6. The survival function, Sss, for the self-similar distribution 

in radius of N/Sl, N/S2, N/S3 and N/S4, as reported in Fig. 4 of 
von Benda-Beckmann fc M iiller (2008) (full red lines), as represented by (|35p . The 
survival function, SF(x, b), of the radius of the distribution which involves the Mei- 
jer G-function for V p (2,3) as represented by (|31[1 when b — 12 Mpc, b — 14 Mpc, 
b — 17 Mpc, and b = 19 Mpc (dashed green lines). 
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TABLE 10 

THE STATISTICAL PARAMETERS OF THE THEORETICAL RADIUS 
OF THE NPVT SPHERES AS REPRESENTED BY EQ. ([3"6] WHEN 
6=31.33 Mpc AND c = 1.768. 



parameter value 

mean 18.23ft. -1 Mpc 

variance 23.31ft. -2 Mpc 2 

standard deviation 4.82ft, -1 Mpc 

skewness 0.072 

kurtosis —0.162 



indicates c = 5, see Eq. ((HJ, or c = 6, the so called Kiang conjecture: we 
will take c as a variable. The resulting distribution in radius once the scaling 
parameter b is introduced is 

4c (4/3^) C " 1 e- 4/3£I ^7ri? 2 
F K (R,b,c) = ^ . (36) 

The average radius is 

-_ ^3&r(l/3 + c) 

H - 2^^r(c) ' {67) 

and the variance is 

_ 3 2/3 2 2/3 6 2 r_ T (2/3 + c) p (c) + ( p (1/3 + c)) 2\ 

a 2 (R) = i /. (38) 

4c 2/3 7r 2/3(p( c ))2 



The skewness is 



(r (c)) 3 c - 3 r (c) r (1/3 + c) r (2/3 + e) + 2 (r (1/3 + c)) 3 



(r(2/3 + c )r( c )-(r(i/3 + c )) 2 ) 



3/2 



(39) 



and the kurtosis is given by a complicated analytical expression. Fig. [7] shows 
a superposition of the effective radii of the voids in SDSS DR7 with a su- 
perposition of the curve of the theoretical PDF in the radius, Fx(R, b, c), as 
represented by Eq. (|3"6"|) . Table [§] shows the values of \ 2 ■ Table [TU] shows the 
theoretical statistical parameters. 

The result of the integration of the fundamental Eq. (fTTj) inserting c=2 
gives the following PDF for the radius of the cuts 

f(r) NPV TK = 3.4148 ^3^Fr2 2 / 3 G^ .V | $$, 1/3i0 ) ( 40 ) 

< r < 1. 
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R (Mpc/h) 

Fig. 7. Histogram (step-diagram) of the effective radius in the SDSS DR7 with a 
superposition of the PDF in radius of the NPVT spheres, Fk(R, b, c), as represented 
by Eq. (36j). The number of bins is 30, b=31.33 Mpc, and c = 1.768. 



TABLE 11 

NPVT PARAMETERS OF f(r) NPVTK, EQ. (|4*Tj) . 



Parameter 


value 


Mean 


0.488 


Variance 


0.0323 


Mode 


0.517 


Skewness 


-.114 


Kurtosis 


2.614 
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R (Mpc/h) 



Fig. 8. Histogram (step-diagram) of the simulated effective radius of SDSS DR7 
with a superposition of the PDF in radius of the PVT spheres as represented by 
Eq. (|3tj[) . The artificial sample has a minimum value of 10/h Mpc, the number of 
bins is 30, 6= 31.5/h Mpc, and c = 1.3. 



The statistics of NPVT cuts with c—2 are shown in Table fTTl 

On introducing the scaling parameter b, the PDF which describes the 
radius of the cut becomes 

16 TT 2 X 6 

~9 



/(*,*>) 



NPVTK 



= 3.4148 \/3^x2 2/3 G 



4.0 
2,4 



b 6 



1/6,1/2 



! .12/3,1/3,0 I " ! 

< r < b. 



The SF of the second NPVT case, SFnpvtk, with the scaling parameter 6, 
is 

SF( x ,b) NPVTK = 1 - 0.2845 tf^Gti ^ | ^;^, 1/3 ,o) ^ (42) 

< r < b. 

A careful exploration of the distribution in effective radius of SDSS DR7 re- 
veals that the detected voids have radius > 10/h Mpc. This observational fact 
demands the generation of random numbers in the distribution in radii of the 
3D cells as given by Eq. (f36| with a minimal value of 10/h Mpc. The artificial 
sam ple is generated t h rough a numerical computation of the inverse func- 



tion 



Brandt fc Gowan ( 19981) and displayed in Fig. [51 the sample's statistics 



are shown in Tabic [T2] 



6. CONCLUSIONS 

PVT Statistics The approach as given by the stereology to the PDF in 
radii of the circles which result from the intersection between a plane and 
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TABLE 12 

THE STATISTICAL PARAMETERS OF THE ARTIFICIALLY 
GENERATED RADIUS WITH A LOWER BOUND OF 10 /h Mpc, c = 1.3 

AND b =31.5/h Mpc. 



parameter 


value 


mean 


18.69ft- 1 Mpc 


variance 


22.74ft- 2 Mpc 2 


standard deviation 


4.76ft" 1 Mpc 


skewness 


0.33 


kurtosis 


-0.623 


maximum value 


31.27ft- 1 Mpc 


minimum value 


10ft,- 1 Mpc 



a randomly disposed spheres of radius R is actually limited to the case of 
mono-disperse sphe res of radius R and to a power law with radius oc Rr a 



Blower et al.l |2002). Here adopting the same type of demonstration we simply 
substitute into formula (|11[) a new distribution for the generalized radii, (i?), 
of PVT. The resulting distribution in radii, (r), of the circles of intersection 
involves the Meijer G-function. A first test on this new PDF for the radii 
was performed on the 2dFGRS catalog and the theoretical V p (2, 3) cells were 
compared with other fitting functions, see Fig. [6j 
NPVT Statistics 

Among the infinite number of 3D seeds which are non-Poissonian, we se- 
lected a distribution in volume which follows a Kiang function as given by Eq. 
(U with c w2. 

A careful comparison with the measured effective radii permits us to say 
that the NPVT case here considered is a good model because it can reproduce 
the 3D average radius and the variance, see Table [TO] The model for the 
effective radius of the voids as given by the Kiang distribution in volumes 
with c variable can also be used to generate an artificial sample of the effective 
radius of the voids, see Fig. [8] and Table IT2l 

Acknowledgements 

I would like to thank the anonymous referee for constructive comments on 
the text and Mario Ferraro for positive discussions on the Voronoi Diagrams. 



REFERENCES 

Abazajian, K. N., Adelman-McCarthy, J. K., Agiieros, M. A., Allam, S. S., Allende 
Prieto, C, & et al. 2009, ApJS , 182, 543 



20 



ZANINETTI 



Barrow, J. D. & Coles, P. 1990, MNRAS , 244, 188 

Bernardeau, F. & van de Weygaert, R. 1996, MNRAS , 279, 693 

Blower, J., Keating, J., Mader, H., & Phillips, J. 2002, Journal of Volcanology and 

Geothermal Research, 120, 1 
Brandt, S. & Gowan, G. 1998, Data Analysis: Statistical and Computational Meth- 
ods for Scientists and Engineers (New York: Springer- Verlag) . 
Bratley, P. and Fox, B. L. 1988, ACM Trans. Math. Softw., 14, 88 
Chiu, S. N., Weygaert, R. V. D., & Stoyan, D. 1996, Adv. in Appl. Prob., 28, 356 
Coles, P. 1991, Nature , 349, 288 

Colless, M., Dalton, G., Maddox, S., & et al. 2001, MNRAS , 328, 1039 

Ebeling, H. & Wiedenmann, G. 1993, Phys. Rev. E , 47, 704 

Elyiv, A., Melnyk, O., & Vavilova, I. 2009, MNRAS , 394, 1409 

Ferenc, J.-S. & Neda, Z. 2007, Phys. A , 385, 518 

Ferraro, M. & Zaninetti, L. 2011, Phys. Rev. E , 84, 041107 

Geller, M. J. & Huchra, J. P. 1989, Science, 246, 897 

Goldwirth, D. S., da Costa, L. N., & van de Weygaert, R. 1995, MNRAS , 275, 1185 
Ho, M.-W., James, L. F., & Lau, J. W. 2007, ArXiv e-prints 
Icke, V. & van de Weygaert , R. 1987, A&A , 184, 16 
Ikeuchi, S. & Turner, E. L. 1991, MNRAS , 250, 519 
Kiang, T. 1966, Z. Astrophys. , 64, 433 

Marinoni, C, Davis, M., Newman, J. A., & Coil, A. L. 2002, ApJ , 580, 122 

Meijer, C. 1936, Nieuw Arch. Wiskd., 18, 10 

— . 1941, Proc. Akad. Wet. Amsterdam, 44, 1062 

Melnyk, O. V., Elyiv, A. A., & Vavilova, I. B. 2006, Kinematika i Fizika Nebesnykh 

Tel, 22, 283 
M0ller, J. 1989, Adv. Appl. Probab., 21, 37 

— . 1994, Lectures on Random Voronoi Tessellations. (Lecture Notes in Statistics. 
87) (New York: Springer- Verlag). 

Okabe, A., Boots, B., & Sugihara, K. 1992, Spatial tessellations. Concepts and Ap- 
plications of Voronoi diagrams (Chichester, NY: Wiley) 

Olver, F., Lozier, D., Boisvert, R., & Clark, C. 2010, NIST Handbook of Mathemat- 
ical Functions. (Cambridge: Cambridge University Press.) 

Pan, D. C, Vogeley, M. S., Hoyle, F., Choi, Y.-Y., & Park, C. 2011, ArXiv e- 
prints: 1103.4156 

Patiri, S. G., Betancort-Rijo, J. E., Prada, F., Klypin, A., & Gottlober, S. 2006, 

MNRAS , 369, 335 
Pierre, M. 1990, A&A , 229, 7 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical 
Recipes in FORTRAN. The Art of Scientific Computing (Cambridge: Cambridge 
University Press) 

Schaap, W. E. & van de Weygaert, R. 2000, A&A , 363, L29 

Sobol, I.M. 1967, U.S.S.R. Comput. Math. Math. Phys., 7, 86 

Subba Rao, M. U. & Szalay, A. S. 1992, ApJ , 391, 483 

van de Weygaert , R. 1991, MNRAS , 249, 159 

van de Weygaert, R. 1991, Ph.D. thesis, University of Leiden 

— . 1994, A&A , 283, 361 

— . 2002, |arXiv:astro-ph/0206427| 

— . 2003, Statistics of Galaxy Clustering - Commentary (Statistical Challenges in 
Astronomy), 156-186 



STATISTICS OF THE VOIDS 



21 



van de Weygaert, R. & Icke, V. 1989, A&A , 213, 1 

van de Weygaert, R. & Schaap, W. 2009, in V. J. Martinez, E. Saar, E. M. Gonzales, 
& M. J. Pons-Borderia, eds, Lecture Notes in Physics Vol. 665 (Berlin: Springer- 
Verlag), pp. 291-311. 

von Benda-Beckmann, A. M. & Miiller, V. 2008, MNRAS , 384, 1189 

Voronoi, G. 1907, J. Reine Angew. Math., 133, 97 

— . 1908, J. Reine Angew. Math., 134, 198 

York, D. G., Adelman, J., Anderson Jr., J. E., Anderson, S. F., Annis, J., Bahcall, 

N. A., Bakken, J. A., & et al. 2000, AJ , 120, 1579 
Zaninetti, L. 1991, A&A , 246, 291 
— . 2006, Chinese J. Astron. Astrophys. , 6, 387 
— . 2010, Serbian Astr. Jour., 181, 19 



Full addresses go here 



A. THE MEIJER G-FUNCTION 



In general the Meijer G-function is defined by the following Mellm-Barncs 
type integral on the complex plane, 
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z~ s ds, (43) 



j—n+l j=m-\-l 



where the contour of integration C is arranged to lie between the poles of 
T(a,i + s) and the poles of T(bj + s). The G-function is defined under the 
following hypothesis. 

• 0<m<q,0<n<p, and p < q — 1; 
. z^0; 

• no pair of bj, j distinct, j = 1,2, ... , to differ by an integer or zero; 

• the parameters o, € C and bj £ C are such that no pole of T(bj + s) , j = 
1, 2, . . . , to coincide with any pole of r(a,; + s),i = 1, 2, . . . , n; 

• a,i — bj 7^ 1, 2, 3, . . . for i = 1, 2, . . . , n and j = 1, 2, . . . , m; and 



• if V = Q, th e n the definition makes sense only f or \z\ < 1 , see iMeijer 
(ll93d[l94lh : lHo et all d2007h : loiver et all (l2010h . 



